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Abstract 

A thermal field, which frequently appears in problems of decoherence, pro- 
vides us with minimal information about the field. We study the interaction 
of the thermal field and a quantum system composed of two qubits and find 
that such a chaotic field with minimal information can nevertheless entan- 
gle the qubits which are prepared initially in a separable state. This simple 
model of a quantum register interacting with a noisy environment allows us 
to understand how memory of the environment affects the state of a quantum 
register. 
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A thermal field is emitted by a source in thermal equilibrium at temperature T. The 
thermal field is the field about which we have minimal information, as we know only the 
mean value of the energy [1]]. It arises frequently in problems involving the coupling of a 
system to its environment which is in thermodynamic equilibrium. Recently, Bose et al. 
showed that entanglement can always arise in the interaction of a single qubit in a pure 
state with an arbitrarily large system in any mixed state and illustrated this using a model 
of the interaction of a two-level atom with a thermal field. Using this model, they studied 
the possibility of entangling a qubit with a large system defined in an infinite dimensional 
Hilbert space. The entanglement between the system and the thermal field reduces the 
system to a mixed state when the field variables are traced over. In this paper we address 
the question 'Is it possible for a thermal field, which is a highly chaotic field, to induce 
entanglement between qubits?' The entanglement of an atom and a field may be a totally 
different problem from the entanglement of two atoms (or qubits) by their mutual interaction 
with a chaotic field. 

We consider a quantum register composed of two two-level atoms interacting with a 
single-mode thermal field. We will investigate the entanglement between the two atoms, 
which are initially separable. This simple interaction model of a quantum register with its 
environment allows us to understand how the memory of the environment affects the state 
of a quantum register. 

When a quantum system of two qubits prepared in p s interacts with an environment 
represented by the density operator pE, the system and environment evolve for a finite time, 
governed by the unitary time evolution operator U(t). The density operator for the system 
and environment at time t is 

m = u(t)(p E ®p s )tf(t). (i) 

After performing a partial trace over environment variables, we find the final density matrix 
p s (t) of the quantum system in the following Kraus representation || 
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p s (t)=Tr E p(t)=Y,K,p s kl (2) 

where the Kraus operators K^s satisfy the property: ^^K^K^ = l s . Unitary evolution of 
the quantum system is a special case in which there is only one non-zero term in the operator 
sum (0). If there are two or more terms, the pure initial state becomes mixed. For the mixed 
thermal environment, not only entanglement but also classical correlation of the system and 
environment can make the system evolve from the initial pure state into a mixed one. The 
mutual information between the system and the environment is non-zero if the information 
for the total system is not the same as the total sum of information for each subsystem. A 
non-zero mutual information is due to classical correlation and/or entanglement ||. Thus 
there being non-zero mutual information between the system and environment is a sufficient 
condition for the initial pure system to have evolved into a mixed state. The action of 
projects the system into a pure state of K^pK^ when the initial state is pure. 

The Jaynes-Cumming model, describing the essential physics of the interaction of a 
single-mode radiation field with a two-level atom, has provided an excellent test ground for 
quantum mechanics In this paper, we are interested in the interaction of two two- 

level atoms with a single-mode thermal field, whose annihilation and creation operators are 
denoted by a and a) . In the context of cavity quantum electrodynamics, the study of the 
interaction with atoms has been used to study for the generation of entanglement between 
two atoms when the two atoms are present simultaneously in the cavity || and when the 
two atoms interact consecutively with the cavity |§ . In these studies, the cavity is normally 
prepared in the vacuum and by the superposition of processes involving depositing and not 
depositing one photon into the cavity, the two atoms can evolve into their entangled state. 
However we are interested in the possibility of entanglement via a chaotic thermal field about 
which we have only minimum information. 

For simplicity, we consider the identical two-level atoms 1 and 2 which are coupled to 
a single-mode thermal field with the same coupling constant 7. The ground and excited 
states for the atom i {i = 1,2) are, respectively, denoted by \g)i and |e)j. The cavity mode 
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is assumed to be resonant with the atomic transition frequency. Under the rotating wave 
approximation, the Hamiltonian in the interaction picture is 



Hi 



h l Yl fat + 



i=l,2 



(3) 



where the atomic transition operators are a~ = \g)i(e\ and af = \e)i(g\. 

The density operator for the combined atom-field system follows a unitary time evolution 
generated by the evolution operator, U(t) = exp(—iHt/h). The evolution operator has been 
extensively studied for a two- and three-level atom Jaynes-Cummings model J7|. For the 
interaction of two two-level atoms with a single-mode field, using Hamiltonian (0) and Taylor 
expansions of sine and cosine functions, we find that the analytical form of the evolution 
operator is given in the atomic basis {|ee), \eg), \ge), \gg)} by 
/ 



U(t) 



2 1 2 a{C - 9)a 



(4) 



2 1 2 a{C - <~))tf + 1 -i-faS -i^aS 

—i'jSd^ \ (cos fit + 1) |(cosfk — 1 

—i-ySd^ |(cos£lt— 1) |(cosf2t+r 
\ 2j*at(C-G)& -ija^S -i-ytfS 2 1 2 a\C - Q)a + 1 

where Cl 2 = GP 1 = 2^ 2 {2a)a + 1) and the time-dependent operators C and S are defined by 



—i'jSa 
—ijSa 
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C — Q cos Clt and S = & 1 sin fit. 
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The thermal radiation field with its mean photon number n is a weighted mixture of 
Fock states and its density operator is represented by pE = J2 n Hn]"^) (n\ where the weight 



function P n is 



1 + n ' 



(6) 



We are interested in the evolution of the quantum system hence the time-dependent den- 
sity operator p s (t) is obtained by tracing over the field variables as in Eq.@: p s (t) = 
J2nm Pn{ m \U(t)\n)p s (Q){n\U\t)\m). We denote the ijth element of the matrix (m\(U)\n) by 
UJj" 1 . Using the orthogonality of the Fock state, p s (t) is obtained in the Kraus representation: 
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where the operators are 



Kl = dtag(U™, U™, U™, U%) + U%(\eg){ge\ + h.c.) 
K n 2 = j2U™- x \ee)(s\ + s/2U^- x \s)(gg\ 
Kl = V2U™ +1 \s)(ee\ + V2U™ +1 \gg)(s\ 
k2 = U™- 2 \ee)(gg\ 

K2 = U2r +2 \gg)(ee\ (8) 

with h.c. denoting the Hermitian conjug ate and \s) = (\eg) + \ge))/y/2. Note that 
determines the process that does not change the mean energy of the field while K% and K% 
determine, respectively, one and two photon absorption processes. Kg and kg, respectively, 
describe atomic transitions through one and two photon emission processes. 

Each Kraus operator projects the atomic state into a pure state if the atomic state 
is initially pure. If the atoms are initially in their excited states the only operation that 
projects the system into an entangled state is kg. For a two qubit system described by 
the density operator p, a measure of entanglement can be defined in terms of the negative 
eigenvalues of the partial transposition M: 

£ = -2j2pr. (9) 

i 

where p^ are the negative eigenvalues of the partial transposition of p. When £ = the two 



qubits are separable ||lT| , |T2| and £ = 1 indicates maximum entanglement between them. By 
substituting the initial condition of the atoms, i.e. p s = p s (0) = \ee)(ee\, into Eq.(|7|), we 
find 

p s {t) = A e \ee)(ee\ + 2A s \s)(s\ + A g \gg)(gg\ (10) 

where A e = E n Pn(U™) 2 , A s = £ P n \U^ +1 \ 2 and A g = £„ P n (f/ 4 T +2 ) 2 . The eigenvalues of 
its partial transposition are then p Q = A s and 2p± = (A e +A g )±[(A e +A g ) 2 — AAeAg+AA 2 ] 1 ^ 2 . 



It is obvious that \i Q and \i+ are always positive. The eigenvalue p- becomes negative if and 
only if A s > J A e A g . In order to analyze this case, let us first consider the possibility of 
entangling two atoms when the field is in a Fock state \£), i.e. P n = 5 n ^. Using the values 
of Ulj m in (f|), it is straightforward to prove that | t/^f 1 2 < |C/ff t/ff +2 | regardless of £ and 
interaction time. It is interesting to note that when two atoms are initially in their excited 
states, they cannot be entangled by simultaneous interaction with a Fock state whose energy 
is definite. In fact, the atoms are not entangled even via the vacuum (n = 0) in this case. 
For their interaction with a thermal field, y / ' A e A g has a lower bound J2 n Pn\UiiU£, m \ which, 
extending the analysis for the Fock state, is found to be larger than the A s . This proves 
that when Fock states can not entangle two atoms by simultaneous interaction, neither can 
their classical mixture, for example a thermal state, can do it. 

It has been shown that a single two-level atom and a thermal field can be entangled due to 
their linear interaction, regardless of the temperature of the field 0. If there are two atoms 
interacting with the thermal field, can we still see the entanglement between the field and 
any single atom? As discussed earlier, by interacting with the thermal field, a pure quantum 
system becomes mixed, so we know that the mutual information between the system and 
thermal field should become non-zero. However, this does not mean that there should be 
atom-field entanglement because the mutual information may include classical correlation as 
well. In order to consider the entanglement between a single atom and the field, we first find 
the single atom-field density operator p a _/(t) by tracing the total density operator p(t) over 
the variables of the other atom. The atom is described as a two-dimensional system while 
the field is in infinite-dimensional space; the density operator p a _f(t) is thus defined in a 
2 x oo dimensional space. An analytical form of entanglement is not available for the system 
in 2 x oo space, so that Bose et al project the atom-field state onto a 2 x 2 subspace and then 
compute the entanglement of formation for each of the outcomes. As this projection can 
be done by local actions, it does not increase the entanglement and the result gives a lower 
bound on the entanglement of the entire atom-field system. We have calculated the average 
of entanglements in this way for projected states, and this is plotted in Fig. 1 for various 
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values of the mean photon number n. We see that the atom and field are always entangled 
for t > despite of the presence of the other atom. After a little algebra we can analytically 
prove this as a lower bound of entanglement occurs in the subspace of |e, 1), |e, 0), \g, 1) and 
\g,0) when Pq |t^2i| 2 (^ii ) 2 > 0- The atom-field entanglement is a key to the decoherence 
process. 

More interestingly, we find in Fig. 2 that atoms are entangled by a thermal field when 
one of the atoms is prepared in its excited state and the other in the ground state. If the 
initial state is \eg), the atomic system is represented by the density operator: 

Ps{t) = Y, p nmZ\eg) + V%\ge)){{fig\V% + (ge\U™) 

n 

+ |t/ 1 T"T|ee)(ee| + \U2 2 n+1 \ 2 \gg}(gg\}. (11) 

One of the partial transposition eigen- 

values may be negative when 4 £ nm P„P m H ^iT" 1 1 2 1 ^T^ 1 ! 2 + U™U™U™ m U™ m ) > 0. 
Substituting [7™ m in (|), the left hand side of the inequality is found to be given by 
(1/4) J2nm PiPnSin 2 fl n t sin 2 Q m t/[(2n + l)(2m + 1)] which is always positive. It is surprising 
that a thermal state, which is a highly chaotic state in an infinite dimensional Hilbert space, 
can entangle two qubits depending on their atomic initial preparation. 

Let us now consider the situation when the atoms are initially both in their ground states. 
In this case, K 2 in Eq.(|8]) is the only operator that projects the atoms into an entangled 
state. The density operator for the atoms has the same form as ( |I~0D but with different 
parameters: A e = E n Pn(U™~ 2 ) 2 , A s = £ P n \U'^ l \ 2 and A g = E n Pn(U^) 2 . The measure 
of entanglement (|S[) is calculated and plotted in Fig. 3, which clearly shows entanglement 
between the two atoms for some interaction times. As we did for the initial state |ee), we 
again first consider the case of Fock-state interaction, in order to provide insights for the 
analysis of the more complicated case of thermal-field interaction. For the interaction of 
atoms with the Fock state \£), one of the eigenvalues for the partial transposition of p s (t) is 
negative when E e g = \U^ l \ 2 - \U^U[ l ± 2 \ > 0. Substituting the values Ufi™ in Eq.(§ into 
the definition of Eg, we find 
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E t = 1 ~ c 

9 2£-l 



;i + c)-i|£(l + c) 



(12) 



where c = cosfi n _it and d = (2£ — 1)/J£(£ — 1). The atoms are not entangled by the 
vacuum interaction as Ei = for £ = 0. We find that the two atoms are entangled when 
c > c s = —l + l/(£+d) for any Fock state (£ > 1). This means that it is possible to 
entangle two atoms by the linear interaction with any Fock state of the cavity field if the 
atoms are initially prepared in their ground states, in sharp contrast to the case of the 
initial preparation of atoms in the excited states. In particular, for the Fock state of |1), the 
atoms are always entangled except when cosf^_ii = 1. The critical parameter c s takes the 
minimum value -1 at £ — 1 and maximizes to around —0.76 at £ = 2 then gradually decreases 
as i gets larger. The maximum value of E is 1/ £ at c = —1 + 1/ £ so that entanglement is 
smaller when £ gets larger. 

We now consider the two atoms interacting with a thermal field. When the atoms are 
initially in \eg), entanglement is induced by a thermal field regardless of its mean photon 
number. Is it still true for the atoms initially in \gg)7 We find the answer is no. The 
thermal state is a weighted mixture of Fock states with the weight Pe given by Eq.(|6|). The 
weight factor Pi is a decreasing function with £. When n is very small, because of the 
vacuum dominance the atoms are not very much entangled. As n gets moderately larger, 
Pi becomes more pronounced. Because |1)(1| gives strong atomic entanglement, the atoms 
are entangled more in this case. However, as n ^> 1, the weight function becomes flat and 
cancelation between the component Fock states washes out the entanglement feature. 

So far we have considered the case when the atoms are in their pure states. Bose et al. 
comment that the purity of the qubit is an important ingredient to see entanglement 
between the qubit and a massive system but did not show the dynamics of the entanglement 
when the qubit is initially prepared in its mixed state. It is not possible to analyze the 
atom-field entanglement when the atoms are in a mixed state because of the lack of tools to 
analyze such a mixed 2 x oo system. However, we can analyze how the initial mixedness of 
atoms affects the entanglement between the atoms. We assume that each atom is initially 



S 



prepared in a thermal mixture so that their initial density operators are 

p, = U^ lfl [X\e) i (e\ + (l-X)\g) i {g\] (13) 

where A depends on the temperature of the atoms. Using the Kraus representation (|7|) for 
the initial condition (|TB|), the evolution of the density operator for the two atoms is found in 
the form: p s (t) = A e \ee)(ee\ + A s \s)(s\ + A a \a)(a\ + A g \gg) (gg\ where \a) = (\eg) + \ge))/y/2 
and A:,s,a,g are time-dependent coefficients. The atoms are then entangled when (A s — A a ) — 



2^J A e A g > 0. The measure of entanglement is plotted for the interaction of atoms with a 
low temperature field with n = 1. It is remarkable to see that with even a very small amount 
of mixture with A = 0.065, atomic entanglement is nearly washed out. 

In conclusion, we have demonstrated the very interesting result that two atoms can 
become entangled through their interaction with a highly chaotic system depending on the 
initial preparation of the atoms. This study provides a degree of analytical understanding 
of the decoherence mechanism for a quantum system composed of a few qubits when the 
reservoir with which they interact retains some memory [|13[ . 

This work was supported by the UK Engineering and Physical Science Research Council 
and by the Korean Ministry of Science and Technology through the Creative Research 
Initiatives Program under Contract No. 00-C-CT-01-C-35. 
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FIGURES 




FIG. 1. Single atom-field entanglement against the interaction time when the pair of atoms 
is initially prepared in the state \ee) for n = 0.1 (solid), 1.0 (dotted) and 10.0 (dashed). The 
entanglement always exits for all times t > 0. 
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FIG. 2. Atom- atom entanglement induced by interaction with a thermal field when the atoms 
are initially prepared in \eg) for n = 0.1 (solid), 1.0 (dotted) and 10.0 (dashed). 
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FIG. 3. Atom-atom entanglement induced by interaction with a thermal field of n = 1 when 
the atoms are initially in a pure state ground state (solid), i.e. A = in Eq.([D|), and in a mixed 
state of A = 0.05 (dotted). In the case of A = 0.065, entanglement is shown by a dashed line at the 
spot indicated by a triangle. 
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